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ABSTRACT 

Aims. We investigate the validity ol the mass segregation indicators commonly used in analysing young stellar clusters. 
Methods. We simulate observations by constructing synthetic seeing-limited images of a 1000 massive clusters (10 4 M Q ) 
with a standard IMF and a King-density distribution function. 

Results. We find that commonly used indicators are highly sensitive to sample incompleteness in observational data 
and that radial completeness determinations do not provide satisfactory corrections, rendering the studies of radial 
properties highly uncertain. On the other hand, we find that, under certain conditions, the global completeness can be 
estimated accurately, allowing for the correction of the global luminosity and mass functions of the cluster. 
Conclusions. We argue that there is currently no observational evidence of mass segregation in young compact clusters 
since there is no robust way to differentiate between true mass segregation and sample incompleteness effects. Caution 
should then be exercised when interpreting results from observations as evidence of mass segregation. 
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1. Introduction 

The issue of mass segregation in globular and open clus- 
ters has been discussed in the literature for over 20 years. 
Historically, the first indicator of mass segregation was sim- 
ply that the brightest, most massive cluster members lay 
closest to the cluster core, whereas the faintest, lower- 
mass members fill the whole extent of the cluster area 
( McNamara fc Sekiguchi|1986[ and references therein) . The 
mass segregation issue has since been complemented with 
specific properties that overall quantify the differences in 
the spatial distribution of high- and low-mass stars. The 
most commonly used for this effect are: (1) the radial de- 
pendence of the mass function (e.g ., Moitinho et al. |l997 



Stolte et al1]2006| |Schweizer|2004| |Gouliermis et al.|| 2004 
or luminosi ty function (e.g., Jones fc Staufl'er|1991||de Grij: 



|et al. 2002), (2) the radi al differences in th e ratio of high- 
to low-mass stars (e.g., Hillenbrand 
mass within some characteristic rad 



1997| ), ( 3) the mean 



1988 Hillenbrand & Hartmann 1998), and 



ms (e.g. 



Sagar et al 
) the m ean 



radius of the two distributions (e.g., Sagar et al. 1988 |de 



Grijs et al. 20021 or the direct comparison of the cumula^ 
tive radial density distribution for the two su bsamples (e.g.. 
McNamara fc Sekiguchi|l986| |Tadros"s|[2005l ) . 



Even though almost all young clusters present one or 
more of these properties, several authors have shown that 
the timescale for dynamical relaxation is typically longer 



than the clusters' age (e.g., Bonnell fc Daviesfl998 ), imply- 
ing that the observed mass segregation should not be dy- 
namic al in origin. Faster phenomena, such as violent relax- 
ation ( Hillenbrand fc Hartmann||1998 Binney & Tremaine 
1987 1 and the fact that the massive stars have shorter reiax- 



ation times ( Hillenbrand fc Hartmann|1998 and references 
therein) still do not seem to be sufficient for explaining the 
profusion of young clusters presenting these properties or 
the implicit degree of mass segregation. The alternative is a 
primordial origin, in which the distribution of massive stars 
in young clusters must reflect the initial con ditions and th e 
processes involved in cluster formation (e.g., Bonnell|2 000). 

Conversely, some authors have indirectly shown that the 
way some indicator s are presented is not stati s tically accu- 



rate. For example, Mai'z Apellaniz & Ubeda (20051 prove 



that a significant bias is introduced when building the mass 
function in equal A log(M/M Q ) bins, as is done in the liter- 
ature. A statistical bias is also potentially introduced when 
studying the radial properties of a cluster by dividing the 
cluster area in fixed-width or constant-area annuli rather 
than equal- number annuli, as each annulus will contain dif- 
ferent numbers of stars changing the statistical significance 
from one annulus to the next. 

In this paper we propose to investigate the validity of a 
few traditional mass segregation indicators using synthetic, 
10 4 M Q -class clusters. We describe the biases that result 
directly from the binning of the data and then explore the 
incompleteness in observed samples and its consequences 
on the mass segregation indicators. 



2. Simulations 

2.1. Synthetic clusters 

We created 1000 synthetic clusters, each containing 2 x 10 4 
stars (total mass of 1.5 x 1 4 M Q ). Each cluster member 



was assigned a mass from a Salpeter (1955) (M > 0.5 M Q ) 
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Fig. 1. Lefi: |Salpeter| (| 1955 | (M> 0.5 M ) and |Kroupa| ( |2001 1 (M< 0.5 M ) mass function used to generate the synthetic 
clusters. Right: JN ormalized King (19621 surface density profile. 



andlKroupal (|2001b (0.01 M < M < 0.5 M ) IMF (Figure 



d log N 
dlogM 



-1.35 0.5 < M/Mq 
-0.3 0.08 < M/Mq < 0.5 
-0.7 0.01 < M/M Q < 0.08. 



(1) 



The radial position (distance to the centre) of each artificial 
star was dra wn randomly and independently of mass from a 
King (19621 radial surface density profile (Figure [T] right): 



E(r) 



So 



(2) 



1 + {r/r c y 

with a core radius r c of 0.2 pc, using the Monte Carlo 
method of the cumulative distribution function. For a clus- 
ter of 2 x 10 4 stars, this yields a central projected density 
Eo around 10 4 pc~ 2 . In this way the synthetic clusters are 
blind to any correlation between position and mass, i.e., 
they are not mass segregated. We have not set any bound- 
ary conditions regarding the physics of star formation or 
dynamical evolution, except for the underlying assumption 
that the probability distribution function from which the 
positions and masses are drawn are a King profile and a 
Kroupa-Salpeter IMF. We have also not set any high mass 
cutoff. 

2.2. Synthetic observations 

To investigate the impact of incompleteness due to crowd- 
ing - a strong limitation in most studies of mass segregation 
in real (massive) clusters - we have used IRAF mkobject 
to build "seeing limited" images of the synthetic clusters. 
The masses were transformed into K-band luminosities us- 



ing the mass-luminosity relation from Ascenso et al. ( 2007 ) 
for a distance of 3 kpc and no interstellar extinction. This 
configuration produced stars up to magnitude 19. Since the 
Monte Carlo algorithm for the positions only generates the 
r polar coordinates, a value for 9 was assigned to each r 
from a uniform random distribution between and 360°. 
Figure [2] shows three of the clusters obtained in this way. 

These images were treated as actually observed clusters, 
in the sense that they were subjected to a source extrac- 
tion algorithm (IRAF daofind), PSF photometry (IRAF 



allstar), and cuts in photometric errors to produce the fi- 
nal samples. These synthetic observations were only sensi- 
tive to stars up to magnitude 16.5, with only ~ 29% of the 
original sources up to this magnitude being detected. The 
1000 catalogues produced in this way are hence incomplete 
sub-samples of the original synthetic clusters. Since they 
are meant to pose as real observations, we will hereafter 
refer to the (incomplete) sub-samples as observed*, always 
maintaining the asterisk to avoid confusion with actual ob- 
servations that are not presented here. The properties of 
the (complete) clusters originally generated by the simula- 
tions will be labeled as "true", since they refer to all the 
stars. 



3. Results 

We tested the most commonly used mass segregation indi- 
cators on synthetic, non-segregated clusters to investigate 
how the way we approach observational data may influence 
our perception of mass segregation in massive clusters. 

For each indicator we investigate: (1) the results ex- 
pected for a non-segregated cluster, (2) the statistical ef- 
fects of binning, and (3) the effects of incompleteness of 
the sample due to crowding. The first item is measured 
directly from the synthetic clusters and averaged over the 
whole set to produce the expected properties of a "perfect 
cluster" . The second concerns the way the quantities are 
represented and how it may affect the analysis. The third 
is measured on the observed* clusters to explore in which 
ways the incompleteness of the observed samples due to 
crowding can mimic the effects of mass segregation. 

We used the full width at half maximum of the stars 
in the simulated observations (5 pixels) as the (arbitrary) 
unit of length. 



3.1. Slope vs radius 

The variation of the mass function (MF) with radius is 
already a traditional diagnosis tool for mass segregation 



flMoitinho et al.||1997| |de Grijs et al.||2002| |Stolte et al 



2002[ jGouliermis i et al.pfMT 



2005 Bica & Bonatto 



Bonatto 



Bica||2005| [Stolte 



2005 [ ). In a mass segregated 
cluster we expect to find an increase in the number of mas- 
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Fig. 2. Seeing limited images of three simulated clusters. The brightness of the sources corresponds to the K-band. 



sive stars with respect to the number of low-mass stars 
toward the centre, which translates into a flattening of the 
mass function or an increase of the high-mass end slope, 
hereafter referred simply as slope or T. Conversely, in a 
non-segregated cluster we expect the slope to be constant 
with radius. 



3.1.1. Binning effects 

This section discusses the MF slope analysis performed on 
the original 2 x 10 4 -star clusters. 

When investigating the radial dependence of the mass 
function we must bin the data twice: first we divide the 
cluster area into concentric annuli and then bin the masses 
of the objects in each annulus to produce the mass function 
for that annulus. In order to keep the statistical significance 
and avoid biases, these bins should be defined such as to 
keep the number of stars per bin constant. Historically, the 
radial bins are defined as fixed-width or fixed-area annuli, 
whereas the mass bins are defined from the histogram of 
log(M/M Q ), as constant Alog(M/M Q ) bins, neither one 
keeping the number of stars per bin constant. Instead, the 
radial bins should be defined as equal-number annuli, and 
the mass function as a histogram with each bin containing 
the sa me number of stars and divided by the resultant bin 
width ( jMaiz Apellaniz fc Ubeda||2005 1 . 

Figure [3] shows the variation of T with radius for the 
1000 synthetic 2xl0 4 -star clusters, calculated with several 
combinations of radial and mass bins. The three panels 
to the left have the masses appropriat ely binned, accord- 



ing to Maiz Apellaniz & Ubeda (20051), whereas the pan- 
els to the right are produced using the more traditional 
mass functions from fixed Alog(M/M ) — 0.2 histograms 
of log (M/Mq). In the two top panels the radii are binned in 
annuli with equal number (100) of stars, the middle panels 
have fixed-width (0.2 pc) radial annuli, and in the bottom 
panels the slope of the MF is measured (cumulatively) in 
circles. 

The profile in panel a) is unbiased since we guarantee 
the same statistical significance in both the radial and the 
mass bins by keeping the number of stars in each bin con- 
stant. As such, we find that T is constant with radius, as 
expected for non-segregated clusters, and equal to the input 
value for the simulations, -1.35. When we change the sta- 
tistical significance of the radial bins by considering fixed- 
width annuli (panel bj) or circles (panel cj) while keeping 



the statistical significance of the mass bins, we still find the 
behaviour expected of non-segregated clusters. 

Conversely, all panels to the right display odd trends 
not reflecting the conditions set for the simulations. Panel 
d) presents a T that is constant with radius, hence not sug- 
gesting mass segregation, but larger than the input value 
of -1.35, illustrating the intrinsic bias i n characterising a 
MF built fro m fixed Alog(M/M ) bins ( |Maiz Apellaniz fc| 
Ubeda|2005 1. In panels e) and f) this bias conspires to pro- 



duce contradictory behaviours: panel e) shows a flattening 
of the MF outward, whereas panel f) shows a flattening of 
the MF inward. The profile in panel f) is what we would 
expect to find in a mass-segregated cluster, although it only 
appears as a consequence of the mass binning. Furthermore, 
as can be seen from the last bin, the overall mass function of 
the cluster as measured in fixed Alog(M/M ) bins comes 
out shallower than Salpetcr, revealing a fundamental un- 
derlying bias in this representation of the mass function, as 
we built the clusters to be Salpeter in the first place. 

This shows that the mass function slope is robust 
against radial binning, only if the mass function itself is 
built in an unbi ased way, namely using the |Maiz Apell amz 



& Ubeda 



( 2005| method to bin the masses. 



3.1.2. Incompleteness effects 

The effects of incompleteness on the radial distribution of 
the mass function slope were tested on the observed* clus- 
ters. The completeness assessment and tentative corrections 
will be addressed in Sect. [4] Figure [4] shows the radial de- 
pendence of r for these clusters using the same binning 
combinations as described above. The light lines correspond 
to the radial profiles for the "true" clusters from Figure [3] 

In all panels, regardless of the binning in radius or mass, 
we find a flattening of the MF toward the centre of the clus- 
ter, a signature typically attributed to mass segregation. 
These profiles are in all similar to those descri bed in the lit- 
erature as indicative of mass segregation (e.g.,|Brandl et al~ 
[T9961 jStolte et al.||2002[ |2006| |Schweizer||2004| |Couliermis 
|et al.|]2004| |Bonatto fc Bica||2005| ). In our case, since we 
know the "true" clusters are not segregated, the finding of 
this signature in the observed* clusters cannot be regarded 
as evidence of mass segregation in the underlying cluster, 
but rather as a consequence of crowding. In the presence of 
incompleteness, the statistical biases arising from binning 
are largely overcome by the fact that the low-mass stars go 
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Fig. 3. Mass function slope (T) as a function of radius for the synthetic clusters. Panel a): T is measured in equal- number 
(100 stars) radial annuli (the radii mark the centre of the annuli) and defining the mass functions such that all the mass 
bins have the same number of stars. Panel b): T is measured in fixed-width (0.2 pc) radial annuli (the radii mark the 
centre of the annuli) and defining the mass function as in panel a). Panel c): T is measured in concentric circles (the radii 
mark the limits of the circles) and defining the mass function as in panel a). Panel d): T is measured in equal- number 
(100 stars) radial annuli and defining the mass functions as fixed Alog(M/M Q ) histograms. Panel e): Y is measured in 
fixed-width (0.2 pc) radial annuli and defining the mass function as in panel d). Panel f): T is measured in concentric 
circles and defining the mass function as in panel d). 



undetected in the cluster core. The mass-binning effects are 
only observed as a flatter mass function in general. 



3.2. Ratio of high- to low-mass stars 

In any given region of a non-segregated cluster, apart from 
fluctuations, there should be the same proportion of high 
and low-mass objects as imposed by the underlying mass 
function. In particular, the ratio of high- to low-mass stars 
should not be dependent on radius. This is indeed what we 
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Fig. 4. Mass function slope (T) as a function of radius for the observed* clusters (dark symbols) when compared to the 
"true" clusters (light symbols). Panel a): T is measured in equal-number (100 stars) radial annuli (the radii mark the 
centre of the annuli) and defining the mass functions such that all the mass bins have the same number of stars. Panel 
b): r is measured in fixed-width (0.2 pc) radial annuli (the radii mark the centre of the annuli) and defining the mass 
function as in panel a). Panel c): T is measured in concentric circles (the radii mark the limits of the circles) and defining 
the mass function as in panel a). Panel d): T is measured in equal-number (100 stars) radial annuli and defining the mass 
functions as fixed Alog(M/M ) = 0.2 histograms. Panel e): T is measured in fixed-width (0.2 pc) radial annuli and 
defining the mass function as in panel d). Panel f): T is measured in concentric circles and defining the mass function as 
in panel d). The profiles for the incomplete, observed* clusters mimic the effects of mass segregation. 



find for the synthetic clusters, regardless of how we divide 
the cluster radially. The light symbols in Figure [5] show this 
profile for a high-mass/low-mass threshold of 10 M Q , and 
radial binning consisting of equal- number annuli (left-hand 
panel), fixed- width annuli (middle panel), and concentric 



circles (right-hand panel). All the profiles are flat, again 
validating the absence of mass segregation in the synthetic 
clusters, and present no signature of statistical biases aris- 
ing from radial binning effects. 
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Fig. 5. Ratio of high- to low-mass stars with radius for a mass threshold of 10 M for the "true" clusters (light symbols) 
and observed* clusters (dark symbols) measured in annuli of fixed number of stars (left), fixed-width annuli (middle) and 
concentric circles (right). The profiles for the incomplete, observed* clusters mimic the effects of mass segregation. 



The dark symbols in the panels show the variation of the 
ratio of high- to low-mass stars in the observed* clusters. 
For all geometries the ratio increases toward the cluster 
core, suggesting an apparent depletion of low-mass stars in 
the core. This is a direct consequence of crowding that does 
not allow for the effective detection of faint sources, rather 
than an actual absence of low-mass stars in the underlying 
cluster that could be imputed to mass segregation. Again, 
this profile is similar to tho se cited in the literature as ev- 
idenc e of mass segregation ( Hillenbrand||1997 Stolte et al. 
20061. 



3.3. Mean mass 

Following the same reasoning as before, the mean mass of a 
non-segregated cluster should be independent of the region 
where we choose to measure it. This is what we find when 
we plot the mean mass in concentric annuli for the synthetic 
clusters (Figure p\ light symbols) , regardless of using fixed- 
number (left-hand panel) or fixed-width (right-hand panel) 
rings. 

Conversely, the observed* clusters (dark symbols) dis- 
play a significant increase of the mean mass toward the 
cluster centre, as the faint stars in the centre are not as 
effectively detected as the massive stars, shifting the mean 
mass to higher values, a signature also often attributed to 
mass segregation. 

3.4. Mean radius of the massive stars 

In the present context we define the mean radius of any 
sample of stars as the mean distance of those stars to the 
centre of the cluster. For each cluster, we measured the 
mean radius of the massive stars and compared it to that 
of the cluster as a whole. The massive star subsample was 
defined using mass thresholds of 1, 5, 10, 15, and 20 M Q . 

We find that both the mean radius of the "true" clusters 
and that of their massive stars have the same value for 
all mass thresholds (Figure [7] left) , although the cluster- 
to-cluster fluctuations increase with increasing threshold. 
This changes for the observed* clusters (right-hand panel), 
as the number of low-mass stars in the centre is significantly 
smaller than before due to crowding, causing the total mean 
radius of the cluster to become larger (1.1 pc), whereas the 
mean radius of the set of massive stars remains roughly 



the same for almost all mass thresholds, as these are not 
affected. 

These profiles match those found in the literature (e.g . , 
Sagar et al.||1988| |Bonnell fc Davies||1998| |Schweizer||2004[ ) , 
where the authors consistently find the massive stars to 
have smaller mean radii when compared to the mean cluster 
radii. 



4. On completeness corrections 

In the previous sections we have shown that incompleteness 
due to crowding will mimic the effects of mass segregation 
in the commonly used indicators. This is so because they 
have the same effect: a depletion of low mass stars in the 
cluster core. The fundamental difference is that, whereas 
mass segregation in young clusters implies a physical pro- 
cess over which the stars of different masses are formed 
or somehow appear spatially segregated, crowding simply 
causes the observer to miss the low-mass stars due to the 
resolution limitation of the instrumental set-up. Many au- 
thors are aware of these limitations and apply more or less 
sophisticated completeness corrections to their samples, but 
how good are these corrections? Since we know in advance 
the exact composition of our synthetic clusters, we used 
one of them to perform a thorough investigation on the 
completeness assessment and correction process. 



4.1. Completeness tests 

The direct comparison of the true and observed* proper- 
ties of a synthetic cluster is the most immediate way to 
gain insight into what is actually lost to observational lim- 
itations. Figure [8] shows the difference between the true 
and observed* brightness of cluster stars as a function of 
distance to the centre, while the colour-code maps the ob- 
served* brightness of the stars. It becomes clear that the 
two relevant consequences of crowding toward the cluster 
core are: (1) hampering source detection due to confusion 
caused by the close proximity of the sources, and (2) in- 
flate the stars' brightness by blending their flux with that 
of unresolved neighbours. As a result, as we move into the 
centre of the cluster, we are less and less sensitive to the 
faint stars, and will tend to overestimate, sometimes by sev- 
eral magnitudes, the brightness of those we do detect. The 
bright stars, on the other hand, are equally detected every- 
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Fig. 6. Mean mass of the stars within annuli of equal number of stars (left) and fixed-width annuli (right) for the "true : 
(light symbols) and observed* (dark symbols) clusters. 
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Fig. 7. Mean radius of the stars with masses higher than the designated mass threshold for the average of the 1000 
"true" clusters (left) and for the observed* clusters (right). The dotted line represents the mean radius of the clusters. 
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Fig. 8. Difference between the true and observed* bright- 
ness for the cluster stars as a function of distance to the 
centre. The colour-code maps the observed* brightness of 
the stars. 



where throughout the cluster and their measured brightness 
is hardly affected by the crowding. 



However insightful as this comparison may be, in real 
clusters we must rely on completeness tests to determine the 
extent to which we may trust the observations, as we lack 
the privileged information about the cluster's true compo- 
sition. To address potential accuracy and reliability issues 
of completeness tests, we computed them for a synthetic 
cluster as if it were an actually observed image. For every 
0.5 magnitude bin, we added artificial stars to the image 
in a grid such that each star is separated from its closest 
neighbour by two times the radius of the PSF +1 pixel. 
By forcing the artificial stars to be in such a grid we sam- 
pled the full extent of the cluster area without adding to 
the crowding. The images for each magnitude were then 
subject to source detection and photometry, and the out- 
put lists of sources were compared with the input grids. 
The completeness for magnitude TOj„ is then defined as the 
number of detected grid stars with measured magnitudes 
m ou t that satisfy the condition \rrii n — m ou t\ < 0.1, divided 
by the total number of stars of magnitude m.; n in the in- 
put grid. The latter condition implies that an artificial star 
blended with a cluster star such that it affects its magnitude 
beyond the reasonable photometric uncertainty is rejected 
for completeness purposes, which happens very frequently 
in the crowded core, mainly for the faint stars. The out- 
come of these tests is therefore a high-fidelity completeness 
assessment that contains information, not only about the 
detection success rate, but also about how blending affects 
the incompleteness. We use the definition of completeness 
described above in the following sections. 
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Fig. 9. Global completeness as a function of magnitude 
from the completeness tests (red solid line), and from direct 
comparison of the observed* and true brightnesses (dotted 
line). 



4.2. Global completeness 

The red solid line in Figure [9] shows the global complete- 
ness - the fraction of artificial stars recovered with respect 
to the input stars in the grid - as a function of magnitude. 
These tests return a 90% global completeness for magni- 
tude 12 (6.2 M Q in our example). The dotted line is the 
true completeness defined here as the fraction of observed* 
cluster stars relative to the true number of stars for each 
magnitude in the synthetic clusteiQ The local disparities 
between the two profiles are due to unresolved (blended) 
sources: whereas blending i s ex cluded for the purpose of 
completeness tests (see Sec. 4.1), it does occur in observa- 
tions - blended stars will appear in the list of observed* 
sources as single stars with good photometry. As a conse- 
quence, blended stars "leak" to different magnitude bins 
and cause the true completeness to be contaminated in an 
unpredictable way. For this reason, and again because the 
completeness tests include only single stars, the profiles fail 
to match for some magnitudes. Nevertheless, the overall 
agreement indicates that the accuracy of the global com- 
pleteness tests is quite reasonable. 



4.2.1. Correcting luminosity functions 

An important and surprising corollary of the validation of 
global completeness tests above is that the global proper- 
ties of the cluster, such as its mass function, can effectively 
be corrected for incompleteness due to crowding. Figure 
10 shows the observed* (solid line), true (dotted line), and 
completeness corrected (red diamonds) luminosity functions 
for this cluster. The latter was derived by dividing the first 
by the completeness profile (red line in Figure [9]), and it is 
indeed very faithful to the true luminosity function for all 
but the last corrected bin, where the correction drops from 




12 14 
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Fig. 10. Comparison of the observed* luminosity function 
(solid line) with the true (dotted line) and completeness 
corrected (red diamonds) luminosity functions. 



23% to 8%. The last reliable bin is four magnitudes fainter 
than the estimated 90% completeness limit (6.2 M Q ), cor- 
responding now to a mass of 0.3 M in our example. This 
implies that one would, in principle, be able to see the first 
break of the mass function in this cluster even though the 
completeness limit is a great deal more massive. 

This example shows the potential of the global com- 
pleteness tests, but we emphasise that it is only valid for 
clusters with the same characteristics as the presented syn- 
thetic clusters, when crowding is the only source of in- 
completeness (e.g., no extinction), and for this particu- 
lar method of evaluating completeness. Most completeness 
studies in the literature, although similar, are not as thor- 
ough as the one described here and must therefore be vali- 
dated before extending this result to other degrees of crowd- 
ing and/or observational configurations. 

4.3. Radial completeness limitations 

The global completeness discussed in the previous section 
describes the average behaviour in the whole cluster area, 
but is not representative of the cluster core where the 
crowdin g is most severe. The completeness tests described 
in Sect. [4?T] can then be analysed in concentric rings about 
the centre of the cluster to estimate the radial dependence 
of completeness. Figure [IT] shows the completeness as a 
function of radius for the different magnitudes. 

Extrapolating from Figure [9] one would trust these 
radial profiles to be a fair representation of the true ra- 
dial completeness. However, when comparing both for any 
given magn itude we instead find a blunt disagreement (see 
Figure 12 for magnitude 14). These differences are en- 



1 A summary of the nomenclature and definitions used here is 
given in Table [Al] of Appendix [A] 



tirely attributable to the blending of unresolved sources: 
when selecting stars of magnitude m from the observed* 
list of sources to assess the true completeness, we include 
(blended) stars that are in reality fainter, while at the same 
time excluding true m th magnitude (blended) stars with in- 
flated (combined) brightnesses. The consequence is that we 
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Fig. 11. Completeness tests as a function of radius for the 
different magnitudes. The line in bold indicates the 90% 
completeness limit determined in Sect. |4.2| 



40 



completeness tests 
true completeness 



40 60 
Radius [FWHM] 



Fig. 12. Completeness as a function of radius from the com- 
pleteness tests (red solid line), and from true completeness 
(dotted line) for magnitude 14. 



appear to be more complete than what the grid complete- 
ness tests suggest simply because we cannot differentiate 
between blended and single sources. 

In terms of completeness corrections, these "magnitude 
leaks" due to blending take much larger proportions than 
they did in the global completeness analysis (Sect 



4.21 



On the global scale the effect of the magnitude leaks from 
the core is largely diluted, allowing for reasonable com- 
pleteness corrections. Conversely, the radial completeness 
tests systematically imply very large (over-) corrections in 
the cluster core, which immediately produce a greatly in- 
flated amount of stars, resulting in our case in an (over- 



corrected cluster with typically 3.5 times more stars (up 
to magnitude 14) than the original synthetic cluster. 

In terms of completeness ass essm ent, both radial com- 
pleteness estimates in Figure 12 agree that: (1) the 
global completeness overestimates the completeness in the 
crowded central regions, and (2) that completeness is 
strongly radially dependent being more severe in the clus- 
ter core, which ultimately confirms the hypothesis that it 
is responsible for the apparent mass segregation. 

To summarise, this analysis shows that there is a radial 
dependence of the completeness affecting primarily the low- 
mass stars that we cannot correct for, so no radial property 
(such as mass segregation) can be legitimately measured in 
the presence of severe crowding. 

5. Conclusions 

We used synthetic non-segregated, compact, and massive 
clusters to investigate the impact of the current approach 
to observational data on mass segregation studies. Our con- 
clusion is that incompleteness due to crowding will produce 
the observed properties of mass segregated clusters, even 
when they are not segregated at all. Crowding causes the 
massive stars to be detected more effectively than the low- 
mass stars, resulting in an apparent depletion of low-mass 
stars in the cluster core, which then produces the char- 
acteristics typically attributed to mass segregation. More 
revealing, radial completeness tests provide erroneous es- 
timates of the incompleteness and, as a consequence, lead 
to severe over-corrections. This is even more unsettling if 
we consider that it is not possible to evaluate the accuracy 
of the completeness determinations with the information 
from the observations alone. This is particularly critical for 
distant, rich clusters or clusters observed with poor spatial 
resolution or sensitivity. 

We have also found that the way to present the data 
may furthermore influence the analysis, although to a much 
lesser extent. In what concerns the radial study of the mass 
function, it is imperative that the slope in each radial an- 
nulus be measured in mas s bins with equal number of stars , 
as described in detail by Mafz Apellaniz & Ubeda (2005). 
If this is so, then the radial binning will not influence the 
analysis. The other indicators (ratio of high- to low-mass 
stars and mean mass of the stars in annuli) are not affected 
by radial binning effects. 

This exercise showed that the study of mass segrega- 
tion cannot be dissociated from an exhaustive and rigorous 
study of completeness - which is not often found in the lit- 
erature - and even then extreme caution must be exercised 
when interpreting radial properties as evidence of mass seg- 
regation. 

The presence of interstellar extinction, not included in 
this analysis, will affect the mass segregation indicators in 
a more unpredictable way. On the one hand, the spatial 
distribution of dust can have all possible geometries, al- 
though it is expected that the massive stars in massive clus- 
ters will clear the dust from the cluster's core much more 
rapidly than they will the peripheries. On the other hand, 
the extinction will affect primordially the fainter, low-mass 
stars, again adding to the incompleteness effect and proba- 
bly contribute, at least in their earlier stages, to aggravate 
the incompleteness problem on the cluster scale. 

An interesting outcome of the completeness analysis is 
that in some cases, such as in the one presented here of 
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a massive cluster with no interstellar extinction, it may 
be reasonable to correct the global luminosity function - 
and, by extension, the mass function - for incompleteness, 
thus gaining information about fainter (low-mass) regimes 
that would otherwise be inaccessible, at least using the 
completeness tests we experimented with. Other tests may 
give different corrections and must be validated beforehand. 
This opens a safe, if not new, door to the study of the break 
of the IMF in massive galactic and extra-galactic clusters. 
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Table A.l. 



Term 



Definition 



true clusters 
observed* clusters 
completeness tests 



true completeness 



Synthetic clusters (Sect. 2.1 1 



Synthetic observations (Sect. 2.2 1. 
Completeness assessment obtained by adding a 
image in a grid. Derived completeness is definec 
magnitudes m out that satisfy the conditio n \rn , 
of magnitude m,„ in the input grid (Sect. 4.1 1. ' 
real data. 

Completeness assessment obtained by direct coi 
completeness is defined as the number of star s i 
number of stars in the true cluster (Sect. 4.2 1. r 
true composition of the cluster. 
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Appendix A: Nomenclature 

The following table summarises the nomenclature used in the text. 



